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We investigate the formation of fragments with strangeness degrees 
of freedom in proton- and heavy-ion-induced reactions at high rel- 
ativistic energies. The model used is a combination of a dynamical 
transport model and a statistical approach of fragment formation. 
We discuss in detail the applicability and limitations of such a 
hybrid model by comparing data on spectator fragmentation at 
relativistic SIS/ GS /-energies. The theoretical results are analyzed 
in terms of spectator fragmentation with strangeness degrees of 
freedom such as the production of single-A — 3 ' 4,5 ffe hypernuclei. 
We provide theoretical estimates on the spectra and on inclusive 
cross sections of light hypernuclei, which could be helpful for future 
experiments on hypernuclear physics at the new GSI- and J-PARC- 
facilities. 
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1 Introduction 



Hypernuclear physics opens the unique possibility to investigate the properties 
of the hyperon-nucleon (YN) and the hyperon-hyperon (YY) interaction (a 
historical overview can be found in Ref. [1]). Information on the strangeness 
sector of the hadronic equation of state is essential for nuclear astrophysics, 
e.g., for the properties of neutron stars [2], and in exploring exotic states of 
finite nuclei, e.g., neutron-rich strange nuclei and exotic di-baryonic systems 
[3]. 
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Hypernuclear production in reactions between heavy nuclei was first theoret- 
ically proposed by Kerman and Weiss [4]. These authors found high energy 
reactions as the best possibility to create exotic finite nuclear systems with 
finite strangeness. Since then this field of research has been extended and 
attracted mainly theoretical interest, see, e.g., Refs. [5-7]. 

The experimental situation has so far been rather poor due to the short life 
time of hypernuclei [8], which impedes the detection of hypernuclei. However, 
the life time of hypernuclei seen in the laboratory is considerably enhanced 
in relativistic collisions due to relativistic effects. Hypernuclear production in 
proton-induced reactions has been experimentally studied also at COSY [9] 
and the comparison with theoretical predictions based on transport equations 
of the Boltzmann type has been successfully performed [10]. More recently, 
concrete experimental proposals with high energy heavy-ion and proton beams 
at GS1 (Darmstadt, Germany) and J-PARC (Japan), respectively, have been 
suggested by T. Saito [11,12]. 

The production of hypernuclei by relativistic protons has been investigated be- 
fore by non-relativistic transport theory [10] and fully quantum mechanically 
for coherent reactions [13]. In this Letter we use a covariant formulation of 
transport theory. The approach is first applied to the production of hypernu- 
clei in relativistic heavy ion collisions, as planned at GSI and FAIR. Secondly, 
we also consider proton-induced reactions, but in the J-PARC energy regime of 
50 GeV. The initial and intermediate non-equilibrium stages are described by 
covariant transport theory based on Quantum-Hadro-Dynamics (QHD) [14]. 
The formation of fragments in the exit channels is treated in the statistical 
multifragmentation model (SMM) [15], including different models of fragment 
formation (apart from evaporation/fission). The SMM has been found to ac- 
count successfully for a large variety of observables in heavy ion induced frag- 
mentation reactions, including mass spectra and momentum distributions. We 
thus continue our previous fragmentation studies [16] by including now as a 
new feature a coalescence scenario for the formation of (light) hypernuclei in 
dynamical transport simulations. 



2 Theoretical description of reactions 



The theoretical description of hadron-nucleus and heavy-ion reactions is based 
on the semiclassical kinetic theory of statistical mechanics [17]. The covariant 
extension of this equation is the relativistic Boltzmann-Uehling-Uhlenbeck 
(RBUU) equation [18] 
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for the f(x,k*) 1-body phase space distribution function for the hadrons. In 
particular, nucleons, hyperons and all resonances up to a mass of 2 GeV, as well 
as mesons, e.g., pions,kaons, are explicitely propagated. In the collision term 
the short-hand notations /j = f(x, k*) for the particle and /j = (1 — f(x, k*)) 
for the hole distributions are used, with the 4- momentum k*^ = {El, k) where 
Ey. = V ' M* 2 + k 2 . The collision integral explicitly exhibits the final state 
Pauli-blo eking while the in-medium scattering amplitude includes the Pauli- 
blocking of intermediate states. The collision term includes 21 mesons and 60 
resonances. We note that all elastic and inelastic processes, such as resonance 
production and absorption with associated meson production, are included. 
In particular, inelastic processes with strangeness production are explicitly 
accounted for, e.g., BB — > BYK, Bn — > YK, BK — > BK (with isospin ex- 
change between K° and K + ), which are important for the production of hy- 
perons. An exact solution of the set of the coupled transport equations for the 
different hadrons is not possible. Therefore, the commonly used test-particle 
method for the Vlasov part is applied, whereas the collision integral is mod- 
eled in a parallel-ensemble Monte- Carlo algorithm. The baryonic mean-field 
is modeled within the non-linear Walecka model (mean-field approximation of 
the QHD) [14]. In the following calculations the NL2 parametrization of the 
non-linear Walecka model has been adopted, which gives reasonable values 
for the compression modulus (K = 210 MeV) and the saturation properties 
of nuclear matter [18]. For the parameters of the collision integral see Ref. 
[19]. The numerical implementation for the solution of equation (1) is that of 
GiBUU [20]. 

The formation of hypernuclei depends on the way the fragmentation process 
between nucleons and, in particular, between nucleons and hyperons is ac- 
counted for in dynamical transport situations. It is well known, that the frag- 
mentation process cannot be dynamically described within transport mod- 
els of Boltzmann-type, since in such approaches only one-body phase-space 
densities are calculated and the dynamical phase-space evolution of physical 
fluctuations is neglected. In situations with long time scales for the fragmen- 
tation process, e.g. spectator dynamics in heavy-ion collision and dynamics in 
hadron-induced reactions, a statistical description of the fragmentation pro- 
cess is applicable, as has been shown in a previous work [16]. 

The statistical multifragmentation model (SMM) [15] has been widely ap- 
plied to multiple fragment production. In its standard version [15], SMM does 
not account for strangeness degrees of freedom. Only recently it has been 
attempted to extend the statistical approach by including hyperons [21], how- 
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ever, this work is still in progress. We thus apply the SMM model only for the 
statistical decay of the excited nucleonic spectators. The formation of specta- 
tor hypernuclei is then modeled through a phase space coalescence between the 
fragments and the hyperons originating in the fireball. If hyperons are inside 
the (fragmenting) spectators and their velocities coincide with those of the 
fragments, they are supposed to form a hyperfragment. This is the basic idea 
of the coalescence model [22] , which has been applied in the past to describing 
hypernuclei, see Refs. [5,6]. In those models one assumes that the hypernuclear 
production cross section factorizes into the cross section for hyperon produc- 
tion and that for the fragment production. Phase space coalescence effects are 
accounted for by a coalescence function, which enters into the hypernucleus 
cross section and can be determined within the density matrix formalism, as 
described in detail in Ref. [5]. From those studies the coalescence velocity R 
can be also extracted. R defines the radius in the velocity space, inside of 
which particles can form a (hyper) cluster. From those studies it turned out 
that the coalescence velocity is less or in the order of the Fermi- velocity, which 
agrees with phenomenological investigations of Ref. [7]. 

In this work, the coalescence velocity R for the formation of hyperfragments 
is a phenomenological parameter, which should, in principle, be adjusted to 
empirical information. However, since there are no precise data on hypernu- 
clear formation available yet, we have chosen R such as to produce results as 
close as possible to predictions of the coalescence models of Refs. [6]. 

We have applied the transport model together with subsequent SMM and 
phase space coalescence calculations to 12 C + 12 C@2 AGeV collisions in spec- 
tator fragmentation and to p + 12 C@50 GeV reactions. The reason for the 
choice of these colliding systems is that they will be experimentally investi- 
gated by the HypHI-collaboration at GSI and JPARC facilities, respectively. 
We thus present here for the first time estimates for such reactions from a 
dynamical model. However, before presenting the results on hypernuclei, a 
benchmark study is first discussed in the next section. 



3 Benchmarks 

According to the coalescence scenario [5] the production cross section of hy- 
pernuclei depends mainly on 3 parameters, the cross sections for strangeness 
and fragment production and the coalescence velocity. In this section we test 
the first two parameters from transport calculations and discuss the time evo- 
lution of spectator properties. 

The application of the combined GiBUU+SMM approach has been found to 
work very well in proton induced reactions in terms of cross sections of global 
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Fig. 1. Time dependence of the average mass number (top-left panel), the excitation 
energy per nucleon (top-right panel), the local baryon density (bottom-left panel) 
and the local pressure (bottom-right panel) of projectile spectators in Au+Au col- 
lisions at 0.4 AGeV incident energy. The different curves at each panel indicate 
calculations at different impact parameters 6 = 6, 8, 10 fm (solid, dashed, dot- 
dashed, respectively). 

mass and charge distributions and, in particular, in terms of neutron num- 
ber distributions of a wide selection of different isotopes [16]. However, in 
heavy-ion collisions (HIC) the situation is more complex with respect to that 
in proton-induced reactions. Depending on the centrality one has to clearly 
separate strongly interacting participant matter ("fireball") from less dynam- 
ically developed one ("spectators"). In a theoretical study this can be easily 
done just from the knowledge of suffered collisions for each particle. However, 
the experimental situation is less trivial, and several phase space limitations 
have to be imposed to separate spectators from fireball sources as precisely 
as possible. In the ALADIN experiment, which we will discuss in this section, 
a rapidity selection has been performed [23]. Particles with y > 0.75y pro j are 
assumed to belong to the spectator sources. In the theoretical calculations we 
have performed the same selection for consistency. In addition a density cut of 
P > Psat/300 (with p sat = 0.16 fm~ 3 being the nuclear saturation density) has 
been applied in order to separate particles bound to spectators from emitted 
ones. 

Fig. 1 shows the time evolution of the projectile-spectator in terms of its mass 
number, excitation energy and pressure and for particular centralities. The 
density p sp and pressure P sp have been extracted at the densest point of the 
spectator region. After its formation, the projectile source stabilizes after the 
freeze-out time of tf r ~ 50 — 70 fm/c as indicated by constant values of mass 
number and excitation energy. Furthermore, at t ~ tf r the local pressure [24] 
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Fig. 2. Impact parameter dependence of the average mass number (left panel) and 
the excitation energy per nucleon (right panel) for projectile-spectator matter. The- 
oretical calculations (black band) are compared with data from the ALADIN col- 
laboration [23] (open diamonds). 

becomes negative and almost constant with time. Also the slope dP sp /dp sp 
becomes negative after t > 55,65, 75 fm/c for b = 6,8, 10 fm, respectively. 
This feature may indicate the onset of spinodal instabilities, assuming a con- 
stant temperature evolution in the center of the spectator. However, fragment 
formation cannot be described in transport theory. We, therefore, proceed 
as follows: The spectator matter reaches local equilibrium at t ~ 40 fm/c, 
corresponding to momentum space isotropy, i.e., equality of longitudinal and 
transverse pressure components. Hence the dynamical evolution has come to 
an end. This slightly excited spectator configuration is then taken as the ini- 
tial state for fragmentation which is treated as a statistical process according 
to the laws of thermodynamics as realized in the SMM approach. The physics 
discussed here is similar also for higher energies (see next section) and has 
been recently applied in describing spallation reactions [16]. 

The calculated parameters of the spectator source are shown in Fig. 2 and 
compared with data from the ALADIN collaboration [23] . The black bands in 
the theoretical calculations indicate the effects of moderate changes on mass 
and excitation energy at the time near the onset of equilibration. The compar- 
ison with the data is reasonable except for the very peripheral events, in which 
the average mass of spectator sources is theoretically underestimated. In this 
respect the ground state description for the initial nuclei within the GiBUU 
approach (or within any other transport model) becomes important. If the ini- 
tial nuclear configuration is not a solution of the Hamiltonian underlying the 
transport model the early stage of the dynamical evolution include artificial 
density oscillations which may lead to spurious particle emission processes. 
Results for very peripheral reactions may be thus affected to some extent. An 
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Fig. 3. Left panel: Relation between the impact parameter and the observable Zb OU nd- 
Right panel: Impact parameter dependence of the mean multiplicity of Interme- 
diate-Mass-Fragments (IMF). Theoretical calculations (black band) are compared 
with data from the ALADIN collaboration [23] (open diamonds). 

improved initialization method in the spirit of a covariant density functional 
approach is under study [25]. 

Applying these parameters of the spectator source to the SMM model, we 
arrive at the results on spectator fragmentation, see Figs. 3. Here the relation 
between the impact parameter and the observable Zb OU nd — J2 z>2 Z (with Z 
being the charge) and the average number of intermediate mass fragments 
(3 < Z < 30) as function of the impact parameter are displayed. The ob- 
servable Zb oun d is well reproduced for all centralities, but less so for the most 
peripheral ones (see discussion above). Interesting is the reasonable description 
of the IMF, even for peripheral collisions. This is an indication of stronger 
evaporation effects in the theoretical calculations; the underestimate in the 
mass distribution is compensated by the slightly higher values of the excita- 
tion energy in the theoretical calculations, see again Fig. 2. 

So far it has been demonstrated that inclusive absolute yields are reproduced 
well in the hybrid GiBUU+SMM approach. However, a correct description 
of velocity distributions of statistically produced fragments is also crucial for 
their subsequent coalescence with hyperons. Fig. 4 shows fragment veloc- 
ity distributions in the projectile frame in comparison with experimental data 
taken from [26] . The velocity distributions are described rather well on a quan- 
titative level, in particular the width of the Gaussian-like fragment velocity 
distributions is well reproduced. 

Similar studies on multifragmentation were performed within nuclear molecu- 
lar dynamics and intranuclear cascade models. More extensive reviews can be 
found in Refs. [27] and [28]. Furthermore, we note that in the past similar hy- 
brid approaches have been applied in dynamical situations. It turned out that 
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Fig. 4. Longitudinal velocity distributions in the projectile frame for 
136 Je + Pb@l AGeV reactions. Theoretical calculations (solid curves) are com- 
pared with experimental data (open diamonds) from [26]. 

momentum distributions are reproduced fairly well in heavy-ion collisions as 
well as in hadron- and pion- induced reactions, see, e.g., Ref. [15], which is an 
important issue before applying coalescence for the formation of hypernuclei. 

For the formation of hypernuclei the production and propagation of strangeness 
degrees of freedom play also an important role. We have tested, that for the 
light C + C-system, which will be studied later, the production cross sec- 
tions of K + mesons are well reproduced (see Fig. 5), although the relativistic 
mean-field model gives a rather strong momentum dependence at these high 
energies. A more detailed discussion of strangeness production can be found 
in Ref. [19,29]. 

We conclude that the combined GiBUU+SMM model describes inclusive mul- 
tifragmentation observables at moderate relativistic energies and also strangeness 
production very well. Since strangeness production is also reproduced it is in- 
teresting to extend this model it to the description of hypernuclear formation 
in high energy heavy-ion collisions and proton-induced reactions. These topics 
will be discussed in the next two sections. 



4 Hypernuclei Formation in heavy-ion collisions 



The production of hypernuclei in energetic collisions between light nuclei is 
one of the major projects being investigated by the HypHI-collaboration at 
GSI. The reason for selecting very light systems is the easier identification of 
hypernuclei via the weak decay of the hyperon into pions. In earlier theoretical 
studies, see, e.g., Refs. [5-7], cross sections of the order of only a few microbarn 
{fib) were predicted, which can be easily understood: in collisions between very 
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E km Cm (GeV) 

Fig. 5. Inclusive invariant production cross section of K + mesons in C + C@2AGeV 
collisions at a polar angle of @i a b = 40°. Theoretical calculations (histogram) are 
compared with data from the KaoS collaboration [30] (open diamonds). 

Table 1 

Inclusive production cross sections for different types of hypernuclei (as indicated) 
for the colliding system 12 C + 12 C@2 AGeV. The contribution from exclusive pion- 
nucleon scattering, i.e. ttN — » YK, is shown separately. 





\H 


{He 


\He 


total yield (fib) 


2.2 


4 


1.4 


pionic contribution (fib) 


0.3 


0.2 


0.03 



light systems, such as 12 C+ 12 C@2 AGeV, secondary re-scattering effects inside 
the spectator matter, important for producing low energetic hyperons, are rare 
processes due to the small interaction volume. The situation is different in 
collisions between heavy nuclei due to the high production rate of strangeness 
and many secondary scattering events. 

The coalescence prescription (in coordinate and momentum space, see again 
section 2) for the formation of hypernuclei has been applied to the description 
of spectator fragmentation in C + C@2AGeV collisions. The separation be- 
tween spectator and participant matter has been adjusted by a rapidity crite- 
rion of \y\ > 0.7y pro jjarg- Inclusive rapidity spectra for different light fragments 
and hyperfragments from spectator matter are shown in Fig. 6. The estimated 
hyperfragment production is ~ 5 — 6 orders of magnitude less than that of 
fragment production in general. This is obvious, since the strangeness produc- 
tion cross sections from exclusive primary channels, i.e. primary BB — > BYK, 
BB — > BBKK, and secondary ones (Bn — > YK and BK — > tcY, B stands for 
a baryon and Y for a hyperon) are very low (orders of few fib) [31]. Among the 
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Fig. 6. Rapidity distributions as function of the rapidity (normalized to the 
projectile rapidity in the cm. frame) of different particle types (as indicated) for 
the system 12 C + 12 C@2 AGeV. 

different processes contributing to the formation of hypernuclei, BB — > BY K 
and KB —>■ AB and the secondary one Bn — > YK give the major contribution 
to the formation of hypernuclei. Strangeness production channels with 3, 4- 
body final states are important in order to produce low energy hyperons that 
can be captured by spectator matter. The same argument also holds for sec- 
ondary re-scattering via elastic hyperon-nucleon and pion-nucleon processes. 



The transport results are summarized in table 1, in which the total inclusive 
hypernuclei production cross sections (first line) and the contributions origi- 
nating from pion-nucleon scattering (inside the spectator matter) are shown. 
The transport calculations predict only moderate contributions to the total hy- 
perfragment cross section from indirect coalescence through the 7r5-channel. 
We have investigated this point in more detail by applying the SMM model to 
different freeze-out times (not shown here) with the result of a rather early for- 
mation of hypernuclei. Thus, we conclude that hypernuclei are mainly formed 
due to the capture of fireball hyperons during the passage stage of the spec- 
tators near the expanding fireball region, with subsequent rescattering of hy- 
perons with spectator matter, i.e. AB — > AB. This results was also confirmed 
by earlier model calculations from Wakai et al. [6], in which at low relativis- 
tic energies around 2 AGeV the coalescence between fragments and hyperons 
from ttB collisions is only of minor importance, which is obvious for collisions 
of light systems at low energies. In principle, antikaons K may also contribute 
to hypernuclear production. However, the formation of hypernuclei through 
antikaons has not been considered here due to their very low production cross 
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Table 2 

Inclusive production cross sections for different types of hypernuclei (as indicated) 
for the system p + 12 C@50 GeV. The contribution from pion-nucleon scattering, i.e. 
ttN — ► YK, is shown separately. 





<?tot \}ib] 


Gpionic 




crtot [lM 


Gpionic 




43 


33 




0.9 


0.8 


\H 


28 


14 


i^e 


35 


16 


\H 


9.5 


2.6 


\He 


9.8 


4.3 


%H 


0.8 


0.2 




16 


1.3 



sections with very high threshold value , e.g. BB — > BBKK. 

We conclude from our exploratory calculations, that hypernucleus formation in 
C + C@2AGeV reactions takes place with estimated production probabilities 
on a few fib. To substantiate this result further theoretical extensions would 
be welcome. For example, a purely statistical treatment of excited fireball-like 
spectators including strangeness content seems possible [21]. 



5 Hypernuclear Formation in high energy p + 12 C reactions 

The colliding system p + 12 C@50GeV will be investigated at the new J- PARC 
facility in Japan [12]. As in the case of low-energy proton induced reactions 
[16] and high energy heavy ion collisions (see previous section), we apply also 
here the same dynamical approach for the pre-equilibrium stage. We note that 
the string fragmentation model in terms of the PYTHIA high energy packages 
[32] is included in the transport simulations. 

As an interesting feature here, transport simulations show a dynamical break- 
up of the excited residual nucleus into a remnant in the target rest frame and 
a moving source. Due to the small interaction region and thus small effect of 
re-scattering inside the initial compound system, the target nucleons which 
collide with the beam protons maintain their high momenta after the colli- 
sion. This finally leads to a pre-equilibrium break up of the initial compound 
system into a moving source, which contains highly excited baryon and also 
anti-baryon states. This feature is shown in Fig. 7 in terms of the temporal 
evolution of the density profiles along the z-axis. This particular situation 
appears even in calculations without using any mean-field potential (cascade 
mode), as one would expect from moderate mean-field effects on the reaction 
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Fig. 7. Density profiles p(0, 0, z) along the beam axis at different time steps as 
indicated. Calculations with (solid) and without mean-field potential (dashed) are 
shown. The considered system is p + 12 C@50GeV at a fixed impact parameter of 
b = 2.5 fm. 

dynamics at so high energies. This is not obvious due to the strong momen- 
tum dependence of the underlying relativistic mean-field model which linearly 
rises with energy, a well known feature of mean-field models of Walecka type 
[33]. In particular, the missing momentum dependence in the cascade calcu- 
lations leads to a moderate density enhancement of the moving source due 
to the missing repulsion of the momentum dependent part of the mean-field. 
However, one should note this uncertainty of the momentum dependence of 
the mean-field at so high energies, in which there is no experimental informa- 
tion available, only moderately affects the yields of produced particles strange 
baryons and pions, in particular the differences between the cascade calcula- 
tions and the full BUU calculations differ only by less than 5%. By considering 
inclusive reactions integrated over the entire centrality region the momentum 
dependent effect decreases. Therefore we do not expect considerable effects on 
the production of hypernuclei due to the uncertainty in the momentum de- 
pendence of the mean-field and continue the discussion on hypernuclei, which 
is the main topic of this work. 



A clear separation between processes relevant for the formation of hypernuclei 
and others is a non-trivial task in reactions at very high relativistic energies. 
However, it is possible to give a rough picture using qualitative arguments. 
Relevant for the formation of hypernuclei are low energy hyperons, as in the 
case of heavy-ion collisions discussed in the previous section. Here processes 
with hyperon production in many-body final states are crucial as well as re- 
scattering of initially produced pions with particles of the moving source. In 
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Fig. 8. Rapidity distribution of different particle types (light fragments and hyper- 
ons) at different time steps, as indicated, for the system p + 12 C@50GeV . 



contrast to the situation of low-energy heavy-ion collisions the produced pion 
yields are now very high, i.e. a few b (compared to the yields in the mb region 
for other particles). 

Concerning the fragmentation mechanism the situation is similar to that of 
the previous section. As described in section 2, non-equilibrium dynamics is 
responsible for a de-excitation in terms of binary processes with multi-meson 
production, e.g., pions, and the achievement of thermal equilibrium for both 
systems, the remnant nuclei at rest and the moving source. We have checked 
in terms of the anisotropy ratio [16] that both systems approach thermal 
equilibrium at t > (15 — 20) fm/c (depending on centrality). Afterwards the 
SMM has been separately applied to fragmentation of both sources, i.e. the 
residual nuclei at rest and the moving source. In the later case the SMM 
model has been applied to moving sources consisting of 4 and more particles, 
otherwise a simple coalescence (in coordinate and momentum space) has been 
adopted. 

The whole situation is illustrated in Fig. 8 in terms of the temporal evolution of 
the rapidity distributions of different particle types (fragments and hyperons), 
as indicated. First, a 2-component structure of the source can be clearly seen 
for the heaviest fragments. The rapidity spectra approach a freeze-out state 
after t > 30 — 40 fm/c, in which the nucleonic content of the moving source is 
enhanced due to, e.g., resonance decays. The time behavior of the A-rapidity 
distribution demonstrates the features discussed above. In particular, some 
of the hyperons are produced inside the moving sources in the same rapidity 
region as that of the fragments. The role of initially produced pions will be also 
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Fig. 9. Rapidity distributions of different particle types for the system 
p + l2 C@50GeV. 

of relevance for the hypernucleus formation (see below). These effects indicate 
a production of hypernuclei inside the moving sources. 

Fig. 9 shows the inclusive rapidity spectra for ife-isotopes, and the corre- 
sponding hypernuclei, respectively. The formation of hypernuclei takes mainly 
place inside the moving sources with energies per nucleon in the range of view 
GeV. Due to extremely high multiplicities of mesons, in particular pions, the 
pion-nucleon scattering contribution to the formation of hypernuclei is now 
the dominant process for the lightest hypernuclei, which it was not the case 
at the lower energies of the previous section. This result, which is again in 
line with earlier theoretical investigations [6], is demonstrated in table 2 for 
the total yields of different hypernuclei (first and third columns) and those 
yields originate only from pion-nucleon scattering. Thus, meson-nucleon re- 
scattering should play an important role in the formation and propagation of 
hypernuclei in highly energetic proton-induced reactions. Other processes of 
similar type may also be of relevance, e.g., KN — > Ytt. 



6 Conclusions and outlook 



We have theoretically investigated the different mechanisms for hypernuclear 
production in reactions relevant to future experiments in hypernuclear physics. 
We have applied an extended GiBUU+SMM transport approach which ac- 
counts for the dynamical pre-equilibrium phase space evolution and also for 
the statistical decay of the asymptotically equilibrated spectator sources lead- 
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ing to spectator fragmentation. The formation of hyperfragments has been 
modeled within a phase space coalescence model. 

The reasonable description of spectator fragmentation at intermediate ener- 
gies and the fragmentation of residual nuclei in proton-induced reactions has 
served as the starting basis to extend the transport studies by accounting for 
fragments with strangeness degrees of freedom (hypernuclei). 

We have applied the model to heavy-ion collisions at intermediate relativistic 
energies and to proton-induced reactions at very high incident energies by 
analyzing the transport calculations in terms of formation of light hypernuclei. 
We have estimated the production cross sections of light hypernuclei in such 
reactions by exploratory transport theoretical calculations. 

This first transport theoretical study on hypernuclear production in reactions 
at low and high incident energies can be further extended. In particular, a fu- 
ture combination of our transport model with the recently developed statisti- 
cal multifragmentation model including strangeness degrees of freedom would 
be a helpful tool. Furthermore, a detailed study on possible influences of the 
hadronic mean-field on hypernuclear formation may be important. Within a 
relativistic framework one can extend the present non-linear Walecka model to 
more general hadronic matter including strangeness degrees of freedom [34], 
which has been partially implemented and is being tested. 

In summary, transport simulations predict significant production cross sec- 
tions for highly energetic hypernuclei. The theoretical estimates may be use- 
ful for the future experiments at GSI and the J-PARC facility, since high en- 
ergy hypernuclei can be easily separated from the background and detected. 
The production of hypernuclei with strangeness S = — 2 may be also of ma- 
jor importance, e.g., S- and double- A- hypernuclei. However, due to very low 
statistics these processes have not been studied in this Letter. We conclude 
that this work provides an appropriate theoretical basis for investigations on 
hypernuclear physics. 
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